Epicyclic frequencies derived from the effective potential: 
simple and practical formulae 
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Abstract. We present and discuss a short and simple derivation of orbital epicyclic 
frequencies for circular geodesic orbits in stationary and axially symmetric space- 
times. Such spacetimes include as special cases analytically known black hole Kerr 
and Schwarzschild spacetimes, as well as the analytic Hartle-Thorne spacetime and 
all numerically constructed spacetimes relevant for rotating neutron stars. Our 
derivation follows directly from energy and angular momentum conservation and 
it uses the concept of the effective potential. It has never been published, except for 
a few special cases, but it has already become a part of the common knowledge in 
the field. 
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1. Introduction 

Properties of congruences of nearly circular geodesic orbits in stationary 
and axially symmetric spacetimes are studied because of their funda- 
mental role in the theory of accretion disks around compact objects 
with strong gravity. The Keplerian angular frequency £Ik, the radial 
epicyclic frequency ujr and the vertical epicyclic frequency ojy are the 
most important characteristics of these orbits. Analytic formulae for 
the three frequencies in Schwarzschild, Kerr and Hartle-Thorne metrics 
have been published many times by several authors (see e.g. Wald, 
1984; Okazaki, Kato & Fukue, 1987; Perez, Silbergleit, Wagoner & 
Lehr, 1997; Nowak & Lehr, 1999; Abramowicz & al., 2003) and are 
well-known. 

Here, we recall our unpublished (Abramowicz & Kluzniak, 2000) 
derivation of the epicyclic frequencies in the general case. This fully 
relativistic derivation is remarkably simple: just four short lines of very 
transparent, easy to check algebra. It uses only invariantly defined 
quantities that have obvious physical meaning, and it closely follows 
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the standard Newtonian derivation that is equally short and based on 
the concept of the effective potential. 1 

The standard Newtonian derivation of epicyclic frequencies is briefly 
recalled in section 2, where we also make a few comments on these 
properties of the orbital frequencies that are not commonly appreciated. 
In section 3, for the completeness of presentation, we review some basic 
facts concerning relativistic formalism that are directly relevant to the 
derivation of the orbital frequencies in Einstein's theory. The reader 
who is interested only in our four-line derivation may skip the next two 
sections and start reading from section 4. 



2. Nearly circular geodesic motion: Newton's theory 

In Newton's theory, circular geodesic orbits (i.e., Keplerian circular 
orbits of free particles) occur at locations r = ro, 9 = 9q, where the 
effective potential U e fi has an extremum with respect to both radial r 
and polar-angle 9 coordinates 2 , 

fdu eS \ fdu cS \ l 2 



V dr J V 09 J ' U1 v ' ' 2r 2 



sin 



The derivatives are taken at a constant angular momentum L. The 
gravitational potential 3>(r, 9) is assumed to be stationary, d^/dt = 0, 
axially symmetric d&/d(j) = 0, and to have the equatorial symmetry 
plane at 9 = ir/2. From the last assumption it follows that one can take 
#o = 7r/2, i.e., that there are circular orbits located in the equatorial 
plane. 

Note, that from the solution L = Lk{t) of equation (1) it follows 
immediately that in the equatorial plane the radial distributions of 
Keplerian angular momentum Lk and Keplerian angular velocity £Ik 
are given by 3 , 



L K = ±{r'-) , n K = ^ = ± ^-) . (2) 

Here, the derivatives are taken at the equatorial plane, 6 = tt/2. 

In the linear regime we can consider separately the case of purely 
radial, Sr = r — r , 59 = 0, and purely vertical, 59 = 9 — 9q, 5r = 0, 



1 The same simple method was used in the text book by Wald (1984) in the 
special case of a static, spherically symmetric metric. The general formulae derived 
by Markovic (2000) are complicated and physically unclear — it is not easy to find 
and correct the few missprints that they contain. 

2 We use the standard spherical coordinates [t, r, 6, <f>]. 

3 Obviously, Keplerian orbits are impossible in the region where {d<b/dr) < 0. 
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epicyclic oscillations about circular orbits in the equatorial plane. Take 
a "perturbed" Keplerian orbit which may be slightly not circular, or 
slightly off equatorial plane, and has the same angular momentum as 
the original circular orbit L = Lq, but with a slightly larger energy 
E = Eq + 8E. Because circular orbits occur at extrema of the effective 
potential, the Taylor expansion of U e s near a circular orbit starts from 
the second-order term, 



SU, 



'd 2 U cH 



cff = 77 



8r 2 , 8U, 



'd 2 U eS 



x Or 2 J ' "~ cS 2 ^ do 2 
A small change in the orbital energy 5E may be therefore expressed as 



be 2 



(3) 




1 



Sr 2 +-(8f) 2 , 8E = - 



1 / d 2 U cS 



dB 2 



56 2 +K 2 (56) 2 . (4) 



The energy of the perturbed orbit is obviously a constant of motion, 
and thus 8E = 0, from which it follows that, 



= Sf 



d 2 U cS 
dr 2 



8r + Sf 



= 89 



'd 2 a 



cff 



dO 2 



89 + r 2 86 



(5) 



We are interested in a non-trivial solution 8r / / 80, for which 
the last two equations take the well-known form of a simple harmonic 
oscillator, 

= uj^8r + Sf, = ^86 + 89 (6) 

with the squared eigenfrequencies being the second derivatives of the 
effective potential with respect to proper geodesic distances in radial R 
and vertical V directions. The radial and vertical directions are defined 
in a coordinate independent way that for the particular case of spherical 
coordinates yields dR = dr, dV = rdO. Thus, one may write, 



d 2 U cS 



d 2 U cS 



UJy = 



d 2 U cS 



1 d 2 U cS 



y dR 2 ) \ dr 2 ) ' ^ dV 2 ) r 2 \ d0 2 

From equation (7), and equations (1), (2) it follows that, 



(7) 



, i 2 - o 2 



dhiL 2 K 
d In r 



LOy 



n 2 K + 



1 d 2 $ 



de 2 



(8) 



Several important deductions could be made from (8). Firstly, it is 
obvious that ojr = Qk everywhere if and only if L 2 K oc r, and this 
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together with equation (2) implies $ oc 1/r. Secondly, uy = VLk if 
(but not only if) the gravitational potential is spherically symmetric. 
Therefore, one concludes that in the special but important case of the 
$ = —GM/r potential, all three characteristic orbital frequencies are 
equal and their squares are positive, 

Wj ? = av 2 = = ^ >o. (9) 

Here M is the mass of the gravity source. The equality of the three 
frequencies means that the orbits in the <& = —GM/r potential are 
periodic and closed (a fact known already to Kepler, who discovered 
that they are ellipses). The positive squares imply stability of both 
radial and vertical oscillations, i.e., the dynamical stability of circular 
orbits in the $ = —GM/r potential. 

The third important deduction from (8) is that in the general case, 
with a potential <&(r,0) / —GM/r, the radial epicyclic oscillations 
around Keplerian circular orbits are stable (i.e., uj£ > 0) if and only if 
the Keplerian angular momentum, |L^(r)|, is an increasing function of 
the radius r of these orbits. This statement, of course, is just a special 
case of the well-known Rayleigh stability criterion. 

The last deduction concerns the stability of the vertical epicyclic 
oscillations. If a non-spherical gravitational potential <&(r, 9) has a pos- 
itive second derivative with respect to 9 at the equatorial symmetry 
plane, the vertical oscillations are stable. Only if the second derivative 
is negative and sufficiently large, these oscillations may be unstable. 

Let us illustrate the above general discussion in terms of a specific 
example. The gravitational potential expansion in terms of spherical 
harmonics (with m = because of the axial symmetry) yields, 

tM) = _™_5WM + ..., p 2 (cos<^ We-i). (io) 

The higher order terms are 0(l/r 5 ) because of the equatorial plane 
symmetry. The quadrupole moment Q is negative for an oblate mass 
distribution, and because rotation typically produces a bulge at the 
equator, Q < is the realistic case to consider. Assuming this, and 
neglecting higher multipoles, we deduce from (10) that, 

L 2 K = GMr--^->0, 
K 2 r 



3GQ 



>0, 



dlnL 2 K \ _ 2Mr 2 + 3Q 
dlnr J ~ 2Mr 2 - 3Q 



(11) 
(12) 
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We see from (11) that a negative quadrupole moment never destabilizes 
vertical oscillations. Equation (12) seems to indicate that a negative 
quadrupole moment destabilizes radial oscillations at small enough 
radii 4 



However, let us write \Q\ = qMR 2 , with R* being the equatorial radius 
of the star and q being a dimensionless parameter. A body with fixed 
M and i?* will have the maximal quadrupole if the whole mass is placed 
in an infinitesimally thin ring of matter located at r = R*, 6 = tt/2. It 
is trivial to calculate the quadrupole moment in this case, 



Q = / r 2 P 2 (cos6)dM (14a) 
Jo 

r , rM i 

= r 2 P 2 (cos0) dM = —RlM. (14b) 

L J r=R t: ,8=n/2 Jq 2 



The above formula proves that the destabilization condition (13) is 
impossible to fulfill outside the star, because it implies r < (\/3/2)R* < 
R*. 

One concludes that according to Newton's theory, the quadrupole 
moment cannot destabilize circular orbits (in the equatorial plane) of 
particles orbiting a body. In consciousness of many astrophysicists this 
conclusion constituteses a "proof" that circular orbits around New- 
tonian bodies, are always stable. Thus, the recent finding of Ams- 
terdamski &, al. (2002) that near very rapidly rotating Newtonian 
Maclaurin spheroids there are unstable orbits (with oj 2 ^ < 0) came 
as a surprise. Kluzniak & al. (2001) found that radial epicyclic oscilla- 
tions are destabilised by octupole and higher moments in the harmonic 
expansion of the potentials. 



3. Stationary, axially symmetric spacetimes 

We assume 5 that spacetimes considered here are stationary and axially 
symmetric, which means that they admit two Killing vectors rf and £*, 
which obey, 

4 If the quadrupole is positive, then for r < tq no circular orbits are possible. 
This provides an example to a situation mentioned in the footnote 3. 

5 We use everywhere the (H ) signature and, occasionally, c = 1 = G units. 
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V (l %) = 0, V (i f fc) = 0, £ fc V fc f/j = 7j fc V fc & • (15) 

Here Vj denotes the covariant derivative, and round brackets denote 
symmetrization with respect to indices they embrace. All formulae that 
we use in this contribution are coordinate independent. However, for 
intuitive illustrations, it is convenient to introduce special coordinates 
that follow the time and axial symmetries and closely resemble the 
Newtonian spherical coordinates [t, r, (ft, 6] used in the previous section. 
In these coordinates 

rf = 5\ e = 5V (16) 
5\ being the Kronecker delta, and the metric takes the form, 



ds 2 = (r/rj) dt 2 + 2(7/0 dtdcft + (ft) deft 2 + g rr dr 2 + g ee d9 2 , (17) 

with a notation convention x l y k gik = {xy). The metric does not depend 
on t and (ft, and we assume that it has an equatorial symmetry "plane" , 
9 = 7r/2. The geodesic circular motion is characterized by the four- 
velocity, u l = dx % jds that obeys, 



V*y = 0, u^A^ + ne), u i u k g ik = \. (18) 

Here f2 is the angular velocity, and A~ 2 = (7777) + 2S7 (r/£) + O, 2 (££) > 0. 
From (15) and (18) it follows that the energy defined by £ = (rju) = ut, 
and the angular momentum defined by C = —(£u) = —u^, are constant 
of motion, because they obey 

u k V k £ = 0, u k V k C = 0. (19) 

Obviously, the specific angular momentum defined by t = C/£, is also 
a constant of motion. The angular velocity O and the specific angular 
momentum £ are related by 

The standard definition of the effective potential is, 

U cS = — In (f - 2ig t <t> + fg^) . (21) 
This, together with the identity 1 = g lk uiu k yields, 

\=£ 2 e - 2U - s , or = In £ - U cS . (22) 
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The three metric components that appear in the effective potential 
formula (21) may be invariantly defined in terms of the Killing vecors rf 
and £*, and expressed by three scalar funtions: <I> being the gravitational 
potential, uj being the angular velocity of frame-dragging, and r being 
the circumferential axial radius, 



tt _ -2* = (££) (9 o „\ 

From these definitions one easily recovers the explicit expressions for 
the three scalars, 



" = f2 = -f|y >0 ' e2 * = (w)+^) >0 ' ( 24 ) 

which show that their Newtonian limits are, as they should be, 
Co — ► 0, r — > rsin^. From these limits, from (21), (23), (24), and from 
£ — > L, £ 2 /r 2 <C 1, one concludes that in the Newtonian limit, the 
relativistic effective potential goes to the Newtonian effective potential, 



W eff = i--ln l-2^-^U[/ eff . (25) 




Because £ — > 1 + E, and E <C 1, one sees that equation (22) has the 
correct Newtonian limit, {0 = ln£ — U c g } — ► { E = U c g }. 



4. The four easy pieces 

We now derive simple and practical formulae for both radial and verti- 
cal orbital epicyclic frequencies in four lines of easy, fully transparent, 
algebra. The dot denotes derivative with respect to the proper time 
s, and x denotes either radial r or polar angle coordinate. We use 
the standard definitions for energy, specific angular momentum and 
effective potential, £ = ut, £ = —u^/ut , U e s = —(1/2) \a.(g tt — 2£g tl ^ + 
g2g<f>if>y They h ave been recalled and explained in the previous section. 
As in the Newtonian derivation described in section 2, we consider 
small oscillations with 5£ = 0, S£ ^ 0, 5£ = 0, that occur either in 
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r or in 9 direction, x(s) — xq = 5x, u x = dx/ds = Sx. The first r 
and 9 derivatives of the effective potential are zero, which corresponds 
to an unperturbed circular orbit at the equatorial plane. The Taylor 
expansion starts from the second term, 5x 2 , and also ends there in the 
lowest order. Our four easy pieces consist of 



1 = u t u t g tt + 2u t u (j> g t ^ + u 4> u 4> g H + u x u x g xx , 

2 l-K^""") (fa)2+ t (ii)2 ' 



-2 



5£_ 



(Sx) 



'^e^Asx^Sx 



(26a) 
(26b) 

(26c) 
(26d) 



The last line, with S£ = 0, Sx ^ 0, has obviously the form of a simple 
harmonic oscilator equation, = lo x 5x + Sx and, with the zeroth order 
version of eq. (26b), this yields the final result, 



dx 2 



(27) 



with dX 2 



~9xxdx 2 > being the proper length in the x direction. 



We see that exactly as in Newton's theory, both radial and vertical 
epicyclic frequencies (squared) are equal to second derivatives of the 
invariantly defined effective potential, with respect to the invariantly 
defined coordinates in radial and vertical directions. 



5. Practical calculations 



We gave a simple derivation of general, physically clear, formulae (27) 
for the two epicyclic oscillation frequencies. In this section we show that 
they are also simple to use in practical calculations, in particular when 
the metric is numerically constructed. Indeed, in order to solve second 
order Einstein's field equations and numerically construct a spacetime 
(e.g., outside a rotating neutron star), one must calculate the metric, its 
first and second derivatives. Thus, one stores during the calculations, 



9 U , g't g H , 9rr, gee, 

n tt t<t> J»t> tt t<t> n 4>4> 

9[r] i 9[ r ] > 9[ r ] > 9 [6] > %] > 9[8] ' 

n tt t<t> n H n tt t<t> n H 

9[rr]' 3[rr]' 9[rr}' 9[e6]' %<?]' 9[ee], 



(28a) 
(28b) 

(28c) 
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where derivatives of a quantity Y are indicated by indices in square 
brackets, dY/dx = Y[ x ], d 2 Y/dx 2 = Yj xx ]. Let us now restrict to the 
equatorial symmetry plane, where all dg lk jdO = g^ = 0. The condition 
for the circular orbit, dU c s/dx = is thus trivially fulfilled for x = 0, 
and for x = r it gives the equation for the Keplerian angular mometum 
for circular orbits, 



a tt Of o'fy _j_ f 2 a i 

9[r] z « 9[ r ] + 1 9[ 







which has the solution, 



y9[ r }) 9[ r ] 9[ r ] 



9 [r] 



(29) 



(30) 



with the plus sign before the square bracket holding for Ik > and 
the minus sign for Ik < 0. 

We use (20), and the radial derivative of eq. (26a) for circular orbits, 
to give two expressions for the orbital Keplerian angular velocity in 
terms of the first derivatives of the metric, 



n 



gtt^K + 9tcj, 



K = ~ 



<p){r]± (9t<p)[ r ] ~ (g<f><t>)[r] (9tt)[ 



9t^K + 



(31) 



Our general result (27), explicitly expressed in terms of the metric 
and its first and second derivatives, takes the form, 



1 4^2^+4^] 

2g rr g tt -2£ K g^ + i 2 K g^ ' 

1 n tt - If 4- P 2 

1 9[00] g^ -l- t K g^ 



LO,. 



+ 



= + 



(32a) 
(32b) 



2g 6 g - 2i K tf* + i\ g<M> 

These frequencies uj r , ujg are measured with respect to the proper time 
of a comoving observer. After dividing by the squared redshift fac- 
tor (u 1 ) 2 = e 2UcS {g u + &k gt4>)~ 2 , one gets the "observed" frequencies 
at "infinity", i.e., those measured with respect to the time t of an 
astronomer at rest very far from the source: 



^ 2 [gtt + ^K9t4>Y 

2g rr 

o2 (gtt + M K g t <t>Y 

\La — - 



2ggg 



9[ rr ] 




+ 4 


y [rr] 


g\ee] 


-M K g$ g] 


+ 4 


U' 
9[ee\_ 



(33b) 
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For an aplication of these formulae to realistic models of rotating neu- 
tron stars, see e.g., Kluzniak & al. (2004). 

By expanding ( gfo - U K + l\ g^ ) [r] = 0, one arrives at 

= {gtt + n K gt<t>) 2 Q 

(-9rr) 

(34) 

Because by construction (30), ±sgn£ K = [(fif^j) 2 - gff g[*]] 1/2 > 0, one 

concludes that f2 2 > if and only if d\£K\/dr > 0. 

Thus, the relativistic Rayleigh criterion for circular Keplerian orbits 
is the same as the Newtonian one [cf. equation (8)]: for stability, the 
absolute value of specific angular momentum should be an increasing 
function of the orbital radius. And, just as in Newtonian theory, the 
epicyclic frequencies can be derived from the second derivatives of an 
effective potential (see eq. [27]). 



t<t>\2 



Y-9l 



' n tt 
9[r] 
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